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The generaUzed Chaplygin gas, which interpolates between a high density relativistic era and a 
non-relativistic matter phase, is a popular dark energy candidate. We consider a generalization of the 
Chaplygin gas model, by assuming the presence of a bulk viscous type dissipative term in the effective 
' thermodynamic pressure of the gas. The dissipative effects are described by using the truncated 

, Israel-Stewart model, with the bulk viscosity coefficient and the relaxation time functions of the 

■ energy density only. The corresponding cosmological dynamics of the bulk viscous Chaplygin gas 

' dominated universe is considered in detail for a flat homogeneous isotropic Friedmann-Robertson- 

M \ Walker geometry. For different values of the model parameters we consider the evolution of the 

cosmological parameters (scale factor, energy density, Hubble function, deceleration parameter and 
luminosity distance, respectively), by using both analytical and numerical methods. In the large time 
limit the model describes an accelerating universe, with the effective negative pressure induced by the 
Chaplygin gas and the bulk viscous pressure driving the acceleration. The theoretical predictions of 
the luminosity distance of our model are compared with the observations of the type la supernovae. 
^ ' The model fits well the recent supernova data. From the fitting we determine both the equation of 

state of the Chaplygin gas, and the parameters characterizing the bulk viscosity. The evolution of 
,1 ^ . the scalar field associated to the viscous Chaplygin fluid is also considered, and the corresponding 

^jq' potential is obtained. Hence the viscous Chaplygin gas model offers an effective dynamical possibility 

for replacing the cosmological constant, and to explain the recent acceleration of the universe. 

PACS numbers: 95.36.+X, 98.80.Es, 04.20.Cv, 95.35.+d 
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CN ; I. INTRODUCTION 

^ ■ 

\ The observations of high redshift supernovae [l| and the Boomerang/Maxima/WMAP data [2], showing that the 
' location of the first acoustic peak in the power spectrum of the microwave background radiation is consistent with the 
^T. , inflationary prediction Vl — 1, have provided compelling evidence for a net equation of state of the cosmic fluid lying in 
^ . the range —1 < w = p/ p < —1/3. To explain these observations, two dark components are invoked: the pressureless 
' cold dark matter (CDM) and the dark energy (DE) with negative pressure. CDM contributes flm 0.3, and is mainly 
: motivated by the theoretical interpretation of the galactic rotation curves and large scale structure formation. DE is 
' assumed to provide Qde ^ 0.7 and is responsible for the acceleration of the distant type la supernovae. There are a 
huge number of candidates for DE in the literature (for recent reviews see [1 and (3). 

One possibility are cosmologies based on a mixture of cold dark matter and quintessence, a slowly-varying, spatially 
inhomogeneous component [5|. An example of implementation of the idea of quintessence is the suggestion that it 
is the energy associated with a scalar field Q with self-interaction potential V{Q). If the potential energy density is 
greater than the kinetic one, then the pressure p — /2 — V{Q) associated to the Q-field is negative. Quintessential 
cosmological models have been intensively investigated in the physical literature Q . 

A different line of thought has been followed in 0, d, [1| , where the conditions under which the dynamics of a self- 
interacting Brans— Dicke (BD) field can account for the accelerated expansion of the Universe have been analyzed. 
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Accelerated expanding solutions can be obtained with a quadratic self-coupling of the BD field and a negative coupling 
constant uj [7|]. 

Dissipative effects, including both bulk and shear viscosity, are supposed to play a very important role in the early 
evolution of the Universe. A cosmic fluid (pressureless and with pressure) obeying a perfect fluid type equation of 
state cannot support the acceleration Q . A solution to this problem, and thus avoiding the necessity of a potential 
for the BD field, is to assume that some dissipative effects of bulk viscous type take place at the cosmological scale 
Q . A combination of a cosmic fluid with bulk dissipative pressure and quintessence matter can drive an accelerated 
expansion phase of the Universe and also solve the coincidence probleni (the observational fact that the energy 
density of cold dark matter and of Q-matter should be comparable today) . The dynamics of a causal bulk viscous 
cosmological fluid filled flat homogeneous Universe in the framework of the BD theory was considered in [l^] ■ The bulk 
viscous pressure term in the matter energy-momentum tensor leads to a non-decelerating evolution of the Universe. 

Neither CDM nor DE have direct laboratory observational or experimental evidence for their existence. Therefore 
it would be important if a unified dark matter - dark energy scenario could be found, in which these two components 
are different manifestations of a single fluid A candidate for such an unification is the so-called generalized 

Chaplygin gas, which is an exotic fluid with the equation of state p = —B/p", where B and n are two parameters 
to be determined, ft was initially suggested in [l^l with n — 1, and then generalized in |15] for the case n ^ 1. The 
Chaplygin gas also appears in the stabilization of branes in Schwarzschild-AdS black hole bulks as a critical theory at 
the horizon fld\ and in the stringy analysis of black holes in three dimensions [13] . The Chaplygin equation of state 
can be derived from Born-Infeld type Lagrangians p^ . [T^ . This simple and elegant model smoothly interpolates 
between a non-relativistic matter phase (p = 0) and a negative-pressure dark energy dominated phase. 

The cosmological implications of the Chaplygin gas model have been intensively investigated in the recent literature 
[20| . The Chaplygin gas cosmological model has been constrained by using different cosmological observations, 
like type la supernovae [2l|, the CMB anisotropy measurements [l^, gravitational lensing surveys Q], the age 
measurement of high redshift objects [131 and the X-ray gas mass fraction of clusters |25|. The obtained results are 
somewhat controversial, with some of them claiming good agreement between the data and the Chaplygin gas model, 
while the rest ruling it as a feasible candidate for dark matter. In particular, the standard Chaplygin gas model 
with n = 1 is ruled out by the data at a 99% level [1^. The exact solutions of the gravitational field equations in 
the generalized Randall-Sundrum model for an anisotropic brane with Bianchi type I geometry, with a generalized 
Chaplygin gas as matter source were obtained in [26| . 

The possibility of constraining Chaplygin dark energy models with current Integrated Sachs Wolfe (ISW) effect 
data was investigated in [131 . In the case of a flat universe the generalized Chaplygin gas models must have an energy 
density such that f^c > 0.55 and an equation of state w < —0.6 at 95% confidence level. The extent to which the 
knowledge of spatial topology may place constraints on the parameters of the generalized Chaplygin gas (GCG) model 
for unification of dark energy and dark matter was studied in f28| . By using both the Poincare dodecahedral and 
binary octahedral spaces as the observable spatial topologies, the current type la supernovae (SNe la) constraints on 
the GCG model parameters were examined. An action formulation for the GCG model was developed in [2^, and 
the most general form for the nonrelativistic GCG action consistent with the equation of state has been derived. The 
thermodynamical properties of dark energy have been investigated in ^30| . For dark energy with constant equation of 
state w > —1 and the generalized Chaplygin gas, the entropy is positive and satisfies the entropy bound. Observational 
constraints on the generalized Chaplygin gas (GCG) model for dark energy from the 9 Hubble parameter data points, 
the 115 SNLS Sne la data and the size of baryonic acoustic oscillation peak at redshift, z = 0.35 were examined in 
(3l| . At a 95.4% confidence level, a combination of the three data sets gives 0.67 < B/pJ^" < 0.83 (where po is the 
present day energy density) and —0.21 < n < 0.42, which is within the allowed parameters ranges of the GCG as a 
candidate of the unified dark matter and dark energy. However, the standard Chaplygin gas model {n — 1) is also 
ruled out by these data at the 99.7% confidence level. A geometrical explanation for the generalized Chaplygin gas 
within the context of brane world theories, where matter fields are confined to the brane by means of the action of a 
confining potential, was considered in [33 |. 

The evolution of the Universe contains a sequence of important dissipative processes, including GUT (Grand Unified 
theory) phase transition, taking place at t ~ 10"^'' s and a temperature of about T « 10^^ K, when gauge bosons 
acquire mass, reheating of the Universe at the end of inflation {t w 10^"^^ s), when the scalar field decays into particles, 
decoupling of neutrinos from the cosmic plasma (i « 1 s, T « 10^'^ K), when the temperature falls below the threshold 
for interactions that keep the neutrinos in thermal contact, nucleosynthesis, decoupling of photons from matter during 
the recombination era (i w 10 s, T w 10^ K), when electrons combine with protons and no longer scatter the photons 
etc. [H. 

The first attempts at creating a theory of relativistic dissipative fluids were those of Eckart [131 and Landau and 
Lifshitz [3^. These theories are now known to be pathological in several respects. Regardless of the choice of the 
equation of state, all equilibrium states in these theories are unstable and in addition signals may be propagated 
through the fluid at velocities exceeding the speed of light. These problems arise due to the first order nature of the 
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theory, that is, it considers only first-order deviations from the equiHbrium leading to parabolic differential equations, 
hence to infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle of causality. 
Conventional theory is thus applicable only to phenomena which are quasi-stationary, i.e. slowly varying on space 
and time scales characterized by mean free path and mean collision time. 

A relativistic second-order theory was found by Israel [36] and developed in [10] and [13, HI] into what is called 
"transient" or "extended" irreversible thermodynamics. In this model deviations from equilibrium (bulk stress, heat 
flow and shear stress) are treated as independent dynamical variables, leading to a total of 14 dynamical fluid variables 
to be determined. For general reviews on causal thermodynamics and its role in relativity see [33| and [39] . Causal bulk 
viscous thermodynamics has been extensively used for describing the dynamics and evolution of the early Universe, 
or in an astrophysical context [40| . 

It is the purpose of this paper to consider the effects of a possible existence of a bulk viscosity of the generalized 
Chaplygin gas on the cosmological dynamics of the Universe. The viscous effects are described by using the truncated 
Israel-Stewart theory p^ . By using the Laplace transformation and the convolution theorem, the second order 
differential equation describing the evolution of the Hubble parameter H is transformed into an integral equation. 
The field equations are solved by means of an iterative scheme. Then the general solutions of the equations are 
obtained in a parametric form in the zero, first, second and mth order approximation, and the relevant cosmological 
parameters (scale factor, energy density, Hubble parameter, deceleration parameter etc.) are obtained. The scalar 
field interpretation of the Chaplygin gas is generalized to take into account the viscosity and dissipative effects. 

In order to compare the predictions of the model with the observational data we have fitted the luminosity distance- 
redshift relation with the latest observational data of the type la supernovae. The model fits well these data. From 
the fitting we determine both the equation of state of the Chaplygin gas, and the parameters characterizing the bulk 
viscosity. Even by taking into account the effect of the bulk viscosity, the n = 1 Chaplygin gas models are ruled out 
by the observations. 

The present paper is organized as follows. The physical model and the basic equations are presented in Section 
II. The evolution equation for the Hubble parameter is studied in Section HI, and the behavior of the cosmological 
parameters is obtained. The observational data have been compared with the theoretical predictions of the model in 
Section IV. In Section V we discuss and conclude our results. In the present paper we use a system of units so that 
87rG = c= 1. 



II. GEOMETRY, FIELD EQUATIONS AND CONSEQUENCES 

Perfect fluids in equilibrium generate no entropy and no frictional type heating, since their dynamics is reversible 
and without dissipation. A perfect fluid model is adequate for the description of many processes in cosmology. 
However, real fluids behave irreversibly^ and some processes in astrophysics and cosmology cannot be understood 
except as irreversible processes. An important irreversible effect is bulk viscosity, which typically arises in mixtures, 
either of different species, as is the case of the radiative fluid, or of the same species, but with different energies, as 
in a Maxwell-Boltzmann gas. Physically, in cosmology we can think of bulk viscosity as an internal friction due to 
the different cooling rates in an expanding gas. The dissipation due to bulk viscosity converts kinetic energy of the 
particles into heat, and thus we expect it to reduce the effective pressure in an expanding fluid [33l[39j. 

For a flat homogeneous Friedmann-Robertson- Walker (FRW) with a line element: 

ds^ = df - a^{t) {dx^ + dy^ + dz^) , (1) 

filled with a bulk viscous cosmological fluid the energy-momentum tensor is given by : 

T^ = [p + p + H) u.u'' - (p + n) ^.f , (2) 

where p is the energy density, p the thermodynamic pressure, H the bulk viscous pressure and Ui the four velocity 
satisfying the condition u^m* = 1. The effect of the bulk viscosity of the cosmological fluid can be considered by 
adding to the usual thermodynamic pressure p the bulk viscous pressure H, and formally substituting the pressure 
terms in the energy-momentum tensor by PcS = p + H. The particle and entropy fluxes are defined according to 
iV* = nu^ and = aN^ — (TH^/2^r)u*, with n is the number density, a is the specific entropy, T > is the 
temperature, ^ is the bulk viscosity coefficient and r > is the relaxation coefhcient for transient bulk viscous effect 
(i.e. the relaxation time). The evolution of the cosmological fluid is subject to the dynamical laws of particle number 
conservation N'^.^ = and Gibb's equation Tda = d{p/n) +pd{l/rL) [sl, H^- In the following we shall also suppose 
that the energy-momentum tensor of the cosmological fluid is conserved, that is T^.^. = 0, where ; denotes the covariant 
derivative with respect to the metric. 
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The gravitational field equations together with the continuity equation T!^^, = imply 

= p, 



2H + = -p-U, 



(3) 
(4) 



p + 3ip + p)H = -3i7n, 



(5) 



where H — a/a is the Hubble parameter. 

For the evolution of the bulk viscous pressure 11 we adopt the truncated evolution equation [ssl [s^ . obtained 
in the simplest way (linear in 11) to satisfy the H-theovem (i.e. for the entropy production to be non-negative, 
S'*j = n^/^r > [i3|-The evolution equation for 11 is given in the framework of the truncated Israel-Stewart theory 

by m 



rn + n = -3£,H, 



(6) 



where ^ the bulk viscosity coefficient and t the relaxation time. The truncated equation is a g ood approximation of 
the full causal transport equations if the condition |nd (a'^r/^T) /dt\ « a^H/T holds [H, [sj. In order to close the 
system of equations dH) and (O we have to give the equation of state for p and specify r and ^. 

We assume that the isotropic pressure p of the cosmological fluid obeys a modified Chapylin gas equation of state 

El, 



B 



(7) 



where < 7 < 1 and 0<n<l. B is a. positive constant. 

When 7 = 1/3 and the comoving volume of the Universe is small {p — > 00), this equation of state corresponds to 
a radiation dominated era. When the density is small, p — > 0, the equation of state corresponds to a cosmological 
fluid with negative pressure (the dark energy) . Generally the modified Chaplygin equation of state corresponds to a 

the matter content is pure dust with p = Q. The 



!/(«+!) 



mixture of ordinary matter and dark energy. For p = {B /j) 

1/2 

speed of sound Vg — {dp/ dp) ' in the Chaplygin gas is given by 



7(1 -Kn) - ^. 

P 



(8) 



For the bulk viscosity coefficient and for the relaxation time of the viscous Chaplygin gas we assume the following 
phenomenological laws 



ap 



(9) 



where 0<7< l,a>0 and s > are constants [42|. Eqs. (O are standard in cosmological models, whereas the 
equation for t is a simple procedure to ensure that the speed of viscous pulses does not exceed the speed of light. 

The truncated Israel-Stewart theory is derived under the assumption that the thermodynamical state of the fluid 
is close to equilibrium, that is the non-equilibrium bulk viscous pressure should be small when compared to the local 
equilibrium pressure |n| << p = jp — B/p^. If this condition is violated then one is effectively assuming that the 
linear theory holds also in the nonlinear regime far from equilibrium. However, for a fluid description of the matter, 
the condition ought to be satisfied. 

To see if a cosmological model accelerates or not it is convenient to introduce the deceleration parameter 



dH' 



dt 



1 



p + 3p + 3U 
2p 



37(71 + 1) 
2n 



1 



(71 + 1)7 



3H 
2^- 



(10) 



The positive sign of the deceleration parameter corresponds to standard decelerating models whereas the negative 
sign indicates accelerated expansion. 

By using the assumptions given by Eqs. Q for the bulk viscosity coefficient and the relaxation time, the evo- 
lution equation for the Hubble parameter H for the viscous dissipative Chaplygin gas dominated flat homogeneous 
cosmological models is obtained from the field equations as 



H 



3(7 + l)i7 + 



nB 
IF 



-H 



-2)1-1 



31- 



-H 



2-2s 



H - 



2a 



-BH 



2-2s-2n 



(7+1) 



2a 



H 



4-2s 



9 . 



(11) 
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III. ITERATIVE SOLUTIONS OF THE EVOLUTION EQUATION 



In order to obtain a simpler form of Eq. (jlip we introduce the dimensionless functions 9 and h by means of the 
definitions 



H = (3^a)T^ ^'^^^ (3^a)^ i, s ^, 

and we denote Ao - (B/3"+i) (3^a)-'''+"^/('-'"^ .s ^ 1/2. 
In these variables Eq. pT|) takes the form 

^ + V2 [(7 + 1) /i + nAo/.-i-2„ ^ ;^2(i-.)l |i _ Ao/i^d— ") + (7 + 1) /j^t^-^) _ /^s ^ ^ 1 
dO'^ L J dt^ 2 

Introducing the new variables h = ^Jy and rj = J ^/ydO, respectively, Eq. (jl3p becomes 
^ + V2 (7 + 1) ^ - 2, + V2 (nAo,— + y-^) 



(12) 



(i?7 



- 2Ao2/^ 



+ 2(7+1)2;^-^ = 0,s^ 



1 



(13) 



(14) 



Therefore for a fixed equation of state and known values of s and n the evolution of the viscous Chaplygin gas 
cosmological models is determined by a single numerical parameter Aq. 

Due to the complicated non-linear character of the evolution equation p^ . it is very difficult to obtain exact 
solutions of this equation in the framework of the truncated Israel-Stewart theory. The cosmological model presented 
above, could be robust if the cosmological solutions of equation p4|) . depicting the causal bulk viscous FRW space- 
time, could be studied for an arbitrary range of values of s, n and 7 in the hope of leading to the possibility of correct 
physical description of a well-determined period in the evolution of our Universe. By using the Laplace transform and 
convolution theorem, the differential equation lfT4|) is equivalent with the following integral equation 



F {ri - x) 



y2(nAo2;-"-i+y^-^) 



drj 



2Ao2/^ 



+ 2(7+1)2;^ 



dx + yo (77) 



where 



F (rj — x) 



1 

26 



(15) 



(16) 



(17) 



^ ^ n+y{0)-y*{0) ^ ^ riO) - n_y{0) 



26 



26 



(18) 



and we denoted 6 — ^/j^+2j + 5/ V2, n± = — {'f + 1) /V2 ± 6 and y*(0) — {dy/dri)^^^, respectively. 

The solution of the integral equation (|15[) can be easily obtained by using the method of successive approximations 
or method of iteration to obtain a solution to any desired accuracy. Taking as an initial approximation the solution 
of the linear part of equation , the general solution of the integral equation (fT5|) can be expressed in the first and 
mth order approximation, m G N , as follows 

yi iv) ^ f F (r; - x) [V2 (nAoyo'""' i^) + i^)) Vo i^) - 2Aoy|"'"" i^) + 2 (7 + 1) vt" (^)] dx + yo (77) , (19) 



Vm (??)=/ F{r]- x) 
Jo 



^^ ("Aoy^'l^' i^) + vi-i i^)) y'm-i - 2Ao2/i4 " (^) + 2 (7 + 1) yL\ (x)] dx+y^_^ {rj) , 



(20) 
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(21) 



We can express the iterative solutions of the gravitational field equations for a bulk-viscous fluid filled FRW Universe 
in the framework of the truncated Israel-Stewart theory for s 7^ 1/2 in the following parametric form (in the following 
equations we write a for the variable of integration in order to distinguish it from the independent variable): 



1 



d<7, a = ape 3 ^ , 



(22) 



P = Povir]) ,p = po 



IV iv) 



(23) 




(24) 



1 dy 



2\/2 2/ (77) dr] 



(25) 



n 



Po 



(7-1-1)2/(77) 



Ao 



(77) 



J_dy 
V2dii 



(26) 



where ao, and £^re arbitrary constants of integration, and we denoted po = 3 (3^a)'^^^^ 

In order to solve the evolution equation iteratively we need to chose the initial conditions for the cosmological 
model. The initial value of the function y (77) can be obtained by fixing the initial value of the Hubble function or, 
equivalently, of the density, by using the equation y (0) = (0) / po = p (0) /po- The initial value of dy/dri\^^^ can 

be obtained by fixing the initial value of the deceleration parameter, so that dy/drjl^^^ = — ^2-^/2/3^ [q (0) -f 1]. In 

this way the mathematical initial conditions are fixed by the physical characteristics of the Universe. 

The behavior of the cosmological parameters of the bulk viscous Chaplygin gas filled homogeneous and isotropic 
universe dust universe, with 7 = 0, and, consequently, p — —B/ p"^, are represented, for some fixed values of n and s 
and for different values of Aq in Figs.[l][5l In Fig. [1] the evolution of the scale factor a is represented as a function of 
the dimensionless time 9. 
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FIG. 1: The scale factor a of the bulk viscous Chaplygin gas filled homogeneous and isotropic dust universe (7 = 0) as a 
function of the dimensionless time Q = a^'^'^^^^'H, for n = 0.1, s = 1/4 and different values of Aq: Aq — 0.01 (solid curve), 
Ao = 0.03 (dotted curve), Ao = 0.05 (dashed curve) and Ao = 0.07 (long dashed curve). 



The time variation of the density of the matter is plotted against the time in Fig. [5) In the expanding universe the 
density is a monotonically decreasing function of the cosmic time. 

The behavior of the bulk viscous pressure 11 of the Chaplygin gas is shown in Fig. [3] The negative bulk viscous 
pressure gives a significant contribution to the total negative pressure of the Chaplygin gas. 

The time variation of the bulk viscosity coefficient of the Chaplygin gas is represented in Fig. 21 Similarly to the 
energy density, the bulk viscosity coefficient is a monotonically decreasing function of time. 



7 




FIG. 2: The dimensionless density function r = p/a^^^^ of the bulk viscous Chaplygin gas fiUed homogeneous and isotropic 
dust universe (7 = 0) as a function of the dimensionless time 9 = a^^^^'^'H, for n = 0.1, s = 1/4 and different values of Aq: 
Ao — 0.01 (solid curve), Aq = 0.03 (dotted curve), Ao — 0.05 (dashed curve) and Ao = 0.07 (long dashed curve). 
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FIG. 3: The dimensionless bulk viscous pressure n = Il/a^'^''^"^''^ of the Chaplygin gas for a homogeneous and isotropic dust 
universe (7 = 0) as a function of the dimensionless time 9 = a^^^^~^''H, for n — 0.1, s — 1/4 and different values of Aq: 
Ao = 0.01 (solid curve), Aq = 0.03 (dotted curve), Ao = 0.05 (dashed curve) and Ao = 0.07 (long dashed curve). 




FIG. 4: The dimensionless bulk viscosity coefficient ( — S^/a of the Chaplygin gas for a homogeneous and isotropic dust universe 
(7 = 0) as a function of the dimensionless time 9 = a^^^-^^^'h, for n = 0.1, s = 1/4 and different values of Ao: Ao = 0.01 (solid 
curve), Ao = 0.03 (dotted curve), Ao = 0.05 (dashed curve) and Ao — 0.07 (long dashed curve). 



The time variation of the deceleration parameter q is represented in Fig. [5l In the limit of the large times 5 < 0, 
showing that the viscous Chaplygin gas filled universe experiences an accelerated cosmological dynamics. For large 
values of 9, after experiencing a super-accelerated phase with q < — 1, the viscous Chaplygin gas filled universe ends 
in a de Sitter regime, with q = —1. 

When the bulk viscosity coefficient ^ is proportional to the square root of the density, ^ ~ p^^^, that is, for s = 1/2, 
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FIG. 5: The deceleration parameter q of the bulk viscous Chaplygin gas filled homogeneous and isotropic dust universe (7 = 0) 
as a function of the dimensionless time 9 = a^/''^~^*'t, for n = 0.1, s = 1/4 and different values of Aq: Aq = 0.01 (solid curve), 
Ao = 0.03 (dotted curve), Ao = 0.05 (dashed curve) and Ao = 0.07 (long dashed curve). 




the transformations introduced in Eqs. (I12p cannot be applied. In this case a set of dimensionless variable is given by 

1 1 

'gn+lX 2(i+n) 



H 



B 



, , 3 /3"+i\ 



(27) 



while the dynamics of the Universe is determined by the parameter x ^ l/\/3ct. In these variables and for s = 1/2 
the evolution equation pTjl takes the form 



^ + [(7 + 1 + X) + nh-^-^-] § - xh'-^- + [(7 + 1) X - 1] = 0, .s = 1. 



The transformations h — ^ and r/ — J ^JydO reduces Eq. to 

+ y2[7 + i + x + "y"'""] ^-2xy-" + 2[(7 + i)x-i]y = o. 



(28) 



(29) 



The general behavior of the viscous Chaplygin gas models with s — 1/2 is qualitatively similar to the case s ^ 1/2. 
Therefore we present only the time evolution of the deceleration parameter q, which is shown in Fig. [S] 
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FIG. 6: The deceleration parameter q of the homogeneous and isotropic dust universe (7 = 0) filled with a viscous Chaplygin 
gas, with bulk viscosity proportional to the square root of the energy density (s = 1/2) as a function of the dimensionless time 
6 — t/sqrta, for n — 0.1 and different values of Aq: Ao = 0.01 (solid curve), Aq — 0.03 (dotted curve), Ao = 0.05 (dashed curve) 
and Ao = 0.07 (long dashed curve). 



In the limit of large times the viscous Chaplygin universe with s = 1/2 ends in a super-accelerated state, with 
g« -1.25. 
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IV. COMPARISON WITH OBSERVATIONAL DATA 



From observational point of view fundamental tests of cosmological models can be performed from the study of the 
propagation in a curved space-time of the light emitted by a source in a distant galaxy (like, for example, a supernova), 
and detected on Earth, on a telescope mirror. The luminosity of the source is defined as L = dEe^n/ dtem, that is, the 
luminosity is the total energy emitted by the source in unit time; the suffix em refers to emission. A telescope detects 
a photon flux F = dEj-f.^/ dtj-f.^/ Am , where the suffix rec refers to reception. The flux is the energy detected on the 
telescope mirror surface Am (assumed to be perpendicular to the incident light beam) per unit time interval j43l |. 

An important observational parameter, the redshift z is defined as 1 + z = a^/a, where oq is the present day value 
of the scale factor, which is usually conventionally taken as 1, oq — 1. From the definition of z we obtain da/dt — 



ao/ (1 + zY 



From the definition of the Hubble function we have H = {l/a){da/dt) = —[1/(1 + z)]{dz/dt), 
(1 + z) H . Due to the cosmological expansion the elementary area changes as a? and the 



dz/dt 

which gives dz/dt — 

frequency to of the light is redshifted during the cosmic evolution so that cx 1/a. Therefore F/ L = [l/Atotj (a/oo)^, 
where Atot represents the proper area of a sphere centered in the light source and containing at the time of reception 
the reception point on its surface. The luminosity distance is defined as jl^ 



AttF 



aoTem (1 + 2) , 



(30) 



where rem is the comoving radius. The comoving coordinate rem can be written in terms of an other radial comoving 
coordinate Xem, so that 



Xem = X (flem) 



da 1 dz 



(31) 



In a fiat (k — 0) FRW geometry we have d^i^z) = (1 + z) Xem- The luminosity distance-redshift relation is given 
by Hi 



dL (z) = (1 + ^) 



dz' 



H{z') 



(32) 



In the case of a bulk viscous Chaplygin gas filled universe the luminosity distance dL{z) can be obtained by 
simultaneously solving the following system of differential equations, with z as independent variable 



- 2(1 + z)H^ + 3i/2 = + :^ _ n, 

dz p" 



(33) 



B 

1 + 7 — 



= -3n, 



(34) 



- a (1 + z) p^-^H^ + n = -Zap'H, 
dz 



(35) 



and 



ddL^z) 1 1 + z 

-dL(z) - 



dz 1 + z ' ' H{z)' 
respectively. In order to simplify this system we introduce a set of dimensionless variables, defined as 

DLiz) 



H{z) ^ Hoh{z), p{z) ^ 3H^r{z),n{z) ^ 3H^T^{z),dL{z) 



(36) 



(37) 



we denote A = B/3^^Hq"^^, and choose a so that 3^^ ^aH^^ ^ — 1. Substitution into Eqs. (PS jl - ipS]) transform these 
equations into the form 



2(1 + z)h'^ + 3/1^ = -377- + — - 37r, 
dz r" 



(38) 
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-(1 



dr 
dz 



(1 + 7)r - 



3r" 



= — Stt, 



(39) 



(1 + z) r' 



dz 



-3r'h, 



(40) 



and 



dPUz) 
dz 



1 + z 



1 + z 

Hz)' 



(41) 



respectively. The initial conditions for the system of Eqs. ([55)) - (PT|) are /i(0) = 1, r(0) = 1, 7r(0) = ttq and Dl{Q) = 
0, respectively. In the equation of state one can take 7 = 0. The physical luminosity distance can be written 
as dL{z;\]n;s) = Hq^Dl(z; \\n; s). Once the dimensionless function Dl{z; X;n; s) is known from the numerical 
integration of the system, the fitting with the observational data will fix the numerical values of the parameters 
A, n, s. 

The function d^ (z) can be measured for distant type la supernovae. The luminosity is evaluated by photometry, 
while the rcdshift is evaluated from the spectroscopic analysis of the host galaxy. Each cosmological model has its 
own prediction for the function (z). Therefore the measured di, (z) data are powerful tests of the cosmological 
models, and the luminosity distance can be used to fit the free parameters of the model by using the observational 
results. 

Riess et al. [44] have recently published a new set of 182 gold supernovae, including new Hubble Space Telescope 
(HST) observations, and recalibrations of the previous measurements. We applied the same tests to the Gold2006 
data set [4^], as described in [i^, by taking 7 = 0. Here n and s were adjusted, and fixed values of A=4.0, 4.5, 5.0, 
5.5, 6.0, 6.5, 7.0, 7.5, 8.0 were taken into account. For any probed values of A, there exists a typical value of n, where 
a wide range of s offers a fit to the supernovae data, as shown in Fig [71 As all panels in Fig. [7] show, the acceptable fits 
occur in several disjunct regions, which are slightly separated by s. The best-fit n value was searched in the marginal 
projections of n (0 < s < 2), instead of finding a global minimum. The varying s caused only slight modifications of 
the fit, which is characterized by the marginal projection of s within the 2-a confidence region. 

The preferred n values and the variation with s for the probed A values are represented in Table HI The first column 
represent the values of A we have used for comparison. The best fits for n are represented in the second column, while 
the best fit values for s as well as the allowed range of the parameter is represented in the third column. In the cases 
A — 5.5 and A = 8.0, respectively, a value of s > 1 is required, and the fit quality increases with increasing s. 



A 


n 


s 


4.0 


0.1 


0.8 (s > 0.2) 


4.5 


0.2 


0.8 (s > 0.3) 


5.0 


0.25 


0.8 {s > 0.5) 


5.5 


0.35 


s > 1 


6 


0.4 


1.7 (s > 0.8) 


6.5 


0.5 


s > 0.7 


7 


0.6 


1 (s > 0.7) 


7.5 


0.65 


poor fit for any s 


8 


0.7 


s > 1 



TABLE I: The preferred n values and the variation with s for the probed A values. 



These results show that the model behaves with some complexity in the prediction of the luminosity distance, as 
the parameters do not behave monotonically (with the exception of n). However, the fitting ellipsoids are well-defined 
(see Fig. ^ , and in general the predictions of the model fit well the supernova data. 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the dynamics of a bulk viscous Chaplygin gas filled flat homogeneous and 
isotropic universe. We have derived and formulated the evolution equations of the system, we have considered their 
behavior by using both analytical and numerical techniques, and we have compared the predictions of our model with 



11 



>i=4.5 1=5.0 




FIG. 7: The fit of the luminosity distance - redshift relation of the bulk viscous Chaplygin gas model to the Gold2006 [i^ l 
supernova data. Different panels show models for different A values, the inner (magenta) contours indicate the 1-a confidence 
regions, the outer (turquoise) contours border the 2-a confidence regions. For all plotted values of A an accurate value of n 
with a good fit can be established, while the fit remains acceptable for a wide range of the parameter s. 
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FIG. 8; The parameter values B corresponding to the best fit values (A, n) represented in a logarithmic scale as a function of 
n. 



the supernova data. The most attractive feature of the Chaplygin gas is that it could explain the main observational 
properties of the Universe without appealing to an effective cosmological constant. Generally, the obtained analytical 
and numerical solutions of the gravitational field equations describes an accelerating universe, with the effective 
negative pressure induced by the Chaplygin gas and the bulk viscous pressure driving the acceleration. 

From the equation of state of the Chaplygin gas with 7 = it follows that for the critical values Pc and pc of 
the pressure and density the parameter Wc = Pel Pc is given by Wc = ~B / p^'^^ — He/ Pc- Evaluating this relation at 
the present time when pc = pco gives B = —WcoPcd^^ ~ n(pco)PcO- The Chaplygin gas behaves like a cosmological 
constant for Wco = — 1, which gives the relation between the constant B and the present day value of the bulk viscous 
pressure as 



B = P, 



n+l 
cO 



I n(p,o)1 




n+l 

) 


I n(p,o)l 


PcO 






PcO 



(42) 



where Hq = 3.24 x lO'^^h s-\ 0.5 < ft. < 1, is the Hubble constant [i,[3|. Since 11 (pco) < 0, the presence of the bulk 
viscous effects can significantly increase the value of B. 

By comparing the model with 7 = to the Gold 2006 supernova data, it turns out that a good agreement with 
these observations can be established for a wide range of the power s g (0.2, 2) which occurs in the phenomenological 
laws ([9]), which characterize the bulk viscosity coefficient. The other viscosity parameter a can be obtained from the 



equation 3^ 
we obtain a 



aH, 



2s-l 



(6.2385 X 10-"s- 



1, and by choosing a value for the Hubble parameter. For h = 0.7 {Hq = 2.268 x 10 ^^s ^) 



0.6 



2.8573 X 1052g3^ 



for the above-established range of the parameter s. As for the 



equation of state of the Chaplygin gas, by taking into account the definition of A, A = B/3"Hq"^^, and for the same 
value of the Hubble parameter, the confrontation with supernova data selects the pairs (n, B) represented on Fig. [S] 
Scalar fields are supposed to play a fundamental role in the evolution of the early universe. The Chaplygin gas 
model can be also described from a field theoretical point of view by introducing a scalar field (j) and a self interacting 
potential U{(j)), with the Lagrangian [ll|, [2l|, [4l| 



(43) 



The energy density and the pressure associated to the scalar field 
are given by 



associated to the bulk viscous Chaplygin gas 



and 



'P B 

--c/((/.) = 7P- — + n. 



(44) 



(45) 



respectively. 

The scalar field and the potential can be obtained from the equations 



(t) - 




ndt. 



(46) 
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and 



Uit) 



B 

1-7 P + — -n 



(47) 



respectively, where ipQ is an arbitrary constant of integration. 

The dependence of the potential U (0) on the scalar field </> is represented in Fig. [HI 



. 54 
. 52 
§ 0.5 

D 

.48 
0.46 




1 1.5 2 2.5 3 3.5 4 4.5 



FIG. 9: The potential U{(j)) of the viscous Chaplygin gas associated scalar field as a function of the scalar field (f> for a dust 
universe (7 = 0), n = 0.1, s = 1/4, and for different values of Ao: Ao = 0.01 (solid curve), Ao = 0.03 (dotted curve), Ao = 0.05 
(dashed curve) and Ao = 0.07 (long dashed curve). 



In conclusion, we have found that the viscous Chaplygin gas model offers a real possibility for replacing the effective 
cosmological constant and to explain the recent acceleration of the universe. 
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